N-Jettiness as a Probe of Nuclear Dynamics 
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We propose the use of A-jettiness (tm), a global event shape, as a probe of nuclear dynamics in 
lepton-nucleus collisions. It characterizes the amount of soft radiation between the jet and nuclear 
beam directions. We give the factorization for the 1-jettiness (ri) distribution for the production of 
a single hard jet (J) in lepton-nucleus collisions: I + A(P) — > J(Pj) + X. Each nuclear target gives 
rise to a unique pattern of radiation, affected by nuclear dynamics, that can be quantified by the 
Ti-distribution. Up to power corrections, the ri-distribution allows for a direct measurement of the 
nuclear PDFs. Additional nuclear-dependent effects will be dominated through power corrections 
of size ~ Ql{A) I '(tiPjt) where Q S (A) is a dynamical nuclear scale and Pjt is the transverse 
momentum of the jet. Such nuclear-dependent effects can be studied through a dedicated program 
to measure n -distributions for a range of nuclei and kinematics. We give numerical results for the 
1-jettiness distribution for the simplest case of a proton target at next-to-leading-log accuracy. 

PACS numbers: 12.38.Bx, 12.38.Cy, 12.39.St, 24.85. +p 



Event shape distributions have played a vital role in 
advancing our understanding of various dynamical as- 
pects of QCD. A new global event shape for exclusive 
iV-jet cross-sections, called A^-jettiness (tjv) [HE], was 
recently introduced to veto additional jets in an inclusive 
manner. The definition of tjv is given by 

rjv = £ mmi {^^}, (1) 

k 

where the set {qi} denote reference four- vectors along the 
beam and the N jet directions. The sum over k runs over 
all the hadronic final state particles of momentum p/~ , and 
the Qi are arbitrary normalization constants on the or- 
der of the hard scale in the process. Different choices of 
the set Qi correspond to different definitions of tn with 
different properties [2]. The region of small tjv corre- 
sponds to events with very little radiation between the 
N jets and the beam directions, with the limit tjv — > 
corresponding to N infinitely narrow and well separated 
jets. 

For processes with nuclei in the initial state, N-jettiness 
can provide a quantitative way to characterize the unique 
pattern of final state radiation thereby probing nuclear 
medium effects. N-jettiness can thus be used as a new 
diagnostic tool for studying various aspects of nuclear 
dynamics. Jet tomography is known to be an important 
tool to diagnose properties of the quark-gluon plasma 
(QGP) produced in relativistic heavy ion collisions [3J. 
The interactions of an energetic jet moving through a 
dense medium change the momentum spectrum of jets 
(jet-quenching) leading to the spectacular phenomenon 
of leading particle or jet suppression observed in rela- 
tivistic heavy ion collisions at both RHIC and the LHC 
[4]. The same interaction also induces additional radi- 
ation to alter the radiation pattern between the beam 
and jet directions and to change the overall jet shape. 



Instead of varying jet shape parameters to observe such 
effects [5], one can alternatively use N-jettiness as a way 
to quantify the radiation pattern, while keeping the jet 
definition unchanged. The combination of jet-quenching 
and N-jettiness measurements can provide a comprehen- 
sive jet tomography to diagnose the properties of a dense 
medium, such as QGP or the medium of ordinary nu- 
clei. Similarly, nuclear dynamics in other processes such 
as nuclear deep inelastic scattering (DIS) can be studied 
through N-jettiness distributions. 

In order to compute [TJ [2] , any standard infrared- 
safe jet algorithm can be employed to obtain A" jets and 
the resulting jet momenta are used as the reference jet 
vectors in the set {qi} in Eq. Q. The only information 
needed from the jet algorithm for the calculation of 
are these reference jet momenta. In the limit tn — > 0, 
where the jets are well separated, different jet algorithms 
will give rise to the same set of reference momenta. The 
jet algorithm dependence is thus power suppressed in the 
region of small r^r. The reference momenta {qi} can also 
be obtained directly from a minimization procedure for 
tjv without the use of a jet algorithm [BJ. 

In the N-jettiness formalism, all final state particles are 
grouped either with one of the beam directions or one of 
the jet directions through the minimization requirement 
in Eq. ([I]) . The i-th jet is then defined to be the set of par- 
ticles grouped in the qi-th direction. Note that the con- 
tribution of particles collinear with one of the reference 
vectors qi or particles with soft momenta is suppressed 
due to the dot products {qi ■ Pk} in Eq. (JlJ. The largest 
contributions arise from particles with sizable momenta 
in a direction distinct from the set {qi}. This allows one 
to isolate events with N narrow jets with only soft radi- 
ation between the jets by restricting to a region of small 
tat. One of the advantages of N-jettiness over other event 
shapes such as Thrust is that it allows one to be exclu- 



2 



sive in the number of final state jets while simultaneously 
being insensitive to the contribution from forward beam 
remnants, which are experimentally difficult to control. 
The forward beam remnant contribution to tjv is sup- 
pressed, since these particles are mostly collinear with 
the beam reference momenta in the set {qi}. Further- 
more, tjv quantifies the shape of radiation between the 
N jets as opposed to variables like Thrust which do not 
distinguish the various jet and beam regions. For more 
details on the N-jettiness formalism we refer the reader 
to Refs. [TJ [5] and references within. 

As a first step, in this paper, we consider a specific 
application of N-jettiness: single jet (N = 1) produc- 
tion in lepton-nucleus collisions. Such an analysis can 
be generalized for multiple jets (N > 1) in a straightfor- 
ward manner. Inclusive production of a single jet with 
a high transverse momentum (Pjt) and rapidity (y) in 
lepton-nucleus collisions, £+A(P) — > J(Pj)+X, is a well- 
defined observable and can be systematically calculated 
in the QCD collinear factorization formalism [7J. We set 
up a factorization formalism, based on the Soft Collinear 
Effective Theory (SCET) [8], that in addition gives the 
1-jettiness distribution for such processes. The value of 
Ti reflects the amount of radiation between the jet and 
nuclear beam directions. By studying the T\ distribution 
for a wide range of nuclei, one can systematically probe 
the effect of the nuclear environment on the observed pat- 
tern of radiation. For larger nuclei, the Ti-distribution is 
expected to be broader with the peak position shifted 
toward larger values of t\. This corresponds to the en- 
hanced hadronic activity between the jet and beam direc- 
tions due to nuclear-medium effects in larger nuclei. Such 
a program can be carried out at the proposed Electron- 
Ion Collider (EIC) where a wide range of nuclear targets 
are planned [?J]. In particular, a measurement of the ratio 
of the cross-sections 

dg(A,n,Pj T , y) 
da{A = l,T U Pj T , y y [ ' 

between a larger nuclear target (A) and the nucleon 
target (A=l) can isolate A-dependent nuclear medium 
effects in the three dimensional configuration space 
(ti, Pj T ,y), allowing for systematic studies of dynamical 
nuclear effects. Many of the nuclear-independent theo- 
retical uncertainties are likely to cancel in the ratio in 
Eq.(§. 

Without requiring the detection of the scattered lep- 
ton, the jet transverse momentum Pjt is the only hard 
scale in this process. The tree-level partonic scattering is 
given by lepton-quark scattering through a virtual pho- 
ton exchange i + q — > £' + q. The t\ distribution is just 
proportional to S(ti), corresponding to the infinitely nar- 
row quark jet, and takes the form 

^3^(0) 

o-o<5(n)y^e^/ q/A (x), (3) 



dydPj T dn 



q 



where <7o 



da 



idP, 



is just the tree-level partonic cross- 



mentum. The sum over q runs over the initial quarks 
and antiquarks, e q is the corresponding fractional charge, 
jq/A is the corresponding nuclear parton distribution 
function (PDF) |10) . <to and initial parton momentum 
fraction x are given by 



47ra e 



Q 3 ev 



t 2 



x = eVP ^ (4) 
1-e-vPjr/Q' ( ) 



respectively, with s, i, and u the usual partonic Mandel- 
stam variables, and Q = ^/s denotes the lepton-nucleon 
center-of-mass (CM) energy. Non-zero values of t\ will 
occur only beyond tree level in perturbation theory. For 
T\ <C Pj T , fixed-order perturbation theory breaks down 
due to the appearance of large logarithms of the form 
of a™ In 2 ™ (Pjt/ti). A systematic resummation of such 
logarithms can be performed [1] in the SCET frame- 
work. Furthermore, non-perturbative effects will modify 
the tree- level result in Eq.([3| and these effects will also 
be incorporated. 

We work in the lepton-nucleus CM frame. At lead- 
ing order in the SCET power counting, the factorization 
formula for the Ti-distribution takes the schematic form 



d 3 a 



dydPj T dri 



H ® B A (g> J (g)S, 



(5) 



where <g> represents a convolution, H is the hard function 
encoding physics of the hard scattering, Ba is the nuclear 
beam function encoding the dynamics of the initial 
state nucleus and the collinear beam radiation, J is the 
jet function for the dynamics of collinear modes along the 
jet direction, and S is a soft function for soft radiation 
throughout the event. We refer the reader to Ref. [TT] 
for more details on the beam function and how it differs 
from the jet function. The invariant mass of the nuclear 
beam jet and the final state jet is characterized by the 
scale \[t\Pjt and the soft radiation has momentum of 
size Ti . The relevant scales in the problem are then given 
by 



Pjt > VnPjT > n, 



(6) 



so that H, Ba, J, and S naturally live at the scales 
~ Pjt, fJ-j ~ VnPjr, and fi s ~ n respectively. For 
perturbative values of /ij, the beam function Ba can be 
matched [TT] onto the usual nuclear PDFs f A as 



B A ~l®fA> 



(7) 



section differential in the jet rapidity and transverse mo- 



where the perturbative coefficient I is independent of the 
nucleus. 

According to Eqs. ([5| and Q, for n ^> Aqcd, the 
Ti-distribution is completely determined in terms of the 
nuclear PDFs and the perturbatively calculable functions 
H, I, and S. For r\ ~ Aqcd, the universal soft function 
S becomes non-perturbative. The soft function S is inde- 
pendent of the nuclear target and can thus be extracted 
from the data on lepton-nucleon scattering and used as 
an input for scattering off heavier nuclei. Furthemore, 
for kinematics where /ij ~ U t\ Pj t ~ ^qcd , the beam, 
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jet and soft functions all become non-perturbative. Such 
effects can be important for n ~ Aq CD / Pj T <C Aqcd- 

The leading power formalism is modified through 
power corrections. The power corrections are character- 
ized through ratios of the scales (ijj, (J-j, fJ>s, an d Q S {A). 
Q S (A) is a dynamical scale which characterizes the size 
of power corrections from multiple scattering in a nuclear 
medium and is a function of the nuclear atomic number. 
For the simplest case of the nucleon, Q 2 S {A = 1) ~ ^qcd- 
For larger nuclei Q 2 (A) ~ A u Aq CD where a denotes the 
power law dependence with atomic number A. If there 
is no color transfer between nucleons, a is expected [12] 
to be ~ 1 /3 corresponding to the length of multiple scat- 
terings in the nucleus. When the scales fin or fij are 
of the same order as Q S (A), the power corrections can- 
not be ignored and a new approach such as the Color 
Glass Condensate approach [13] is needed. Q S (A) is also 
sometimes referred as the saturation scale [13] . 

Power corrections to the hard function correspond to 
multiple hard scatterings at the scale /j,h ~ Pj T [14] , 
Such effects generate power suppressed operators which 
give rise to higher twist nuclear beam functions which are 
then matched onto higher twist nuclear matrix elements 
[12]. Thus, these effects can also give rise to nuclear de- 
pendent effects but will be suppressed by the hard scale 
fin ~ Pj T - The dominant nuclear-dependent power cor- 
rections arise from corrections to the OPE in Eq.Q and 
scale as Q S (A)/ /x 2 . These effects correspond to multiple 
scatterings along the nuclear beam at the jet scale with 
typical size 



that ;/-, = ii- 



2, and n , 



-nj, with the 
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A 2 A c 
QCD A 

nP. 



JT 



(8) 



and can be enhanced by the large number of gluons [13] 
in a large nucleus. Precision measurements of the T\- 
distribution for different values of A and Pj T can allow 
for an independent extraction of a and test the expec- 
tation [9] [12] of a ~ 1/3. Power corrections to the jet 
and soft functions will be universal and independent of 
the nuclear target. However, the convolutions of these 
universal power corrections with the nuclear beam func- 
tions can give rise to differences in the t\ -distribution for 
different nuclei. 

By comparing data with the leading-power factor- 
ization formalism in Eqs. ([5| and ([7|, any measured 
non-trivial nuclear modification to the Ti-distribution is 
strong evidence of coherent multiple scattering, and can 
be of interest in the study of multi-parton correlation, 
energy loss in cold nuclear matter, etc. 

We now give some details on the SCET framework and 
the resulting factorization formula. In the lepton-nucleus 
CM frame, the collinear degrees of freedom along the 
nuclear beam have momenta proportional to the light- 
cone vector n A = (1, 0, 0, 1). We define a conjugate light- 
cone vector n A so that n 2 A = n 2 A — 0, tla • tla — 2, and 
Ua = —nA- Similarly, the light-cone four-momentum 
vectors rij and fij associated with the jet direction so 



i>j = 0, nj ■ nj 
jet direction along fij. 

The collinear modes along the nuclear beam and jet 
directions have momentum scalings 



(n A -p,n A -P,V^ A ) ~ P/t(A 2 ,1,A), 
(nj-p,nj-p,p^)~P JT (\ 2 ,l,X), 



(9) 



respectively where A ~ \Jt\I Pjt and p ±A and p ±J aie 
momentum components perpendicular to the beam and 
jet axis respectively. The soft radiation has momentum 
scaling 

{n A ■ p,n A -P,P ±A ) ~ Pjt(A 2 , A 2 , A 2 ). (10) 

The factorization formula for the differential cross-section 
derived using SCET is given by 



d 3 a 



dydP JT dr_ 



- = Co e 2 J dx J ^ J ds.j J dt a 
6 [a 



xH(xQPj T e v ,pij-.fi H ) 
e y Pj T 



Q-e-vPj T 



J q {sj^j) 



-xl qi (-,*o,Mj) fi/A&Hj) (11) 



xS n 



t„ 



sj 



which is the detailed version of the schematic formulas in 
Eqs. ([5| and §f§ and Q a — xQ, Qj = n,j ■ Pj. In deriv- 
ing the above formula we also used the reference vectors 
q A = xPt 1 , where Pa is the initial state nucleus momen- 
tum and qj = (Pj T coshy, Pj T , Pj T sinhy) which is just 
the momentum of a massless jet with transverse momen- 
tum and rapidity Pj T and y respectively. The last two 
arguments of the hard function H denote that renormal- 
ization group (RG) evolution between the scales hh and 
fij is included. Similarly, the soft function is evolved be- 
tween the scales /is and The jet function J q , the 
beam function to PDF matching coefficient I ql , and the 
nuclear PDF fi/A are all evaluated at the common jet 
scale fij. The tree- level values for various functions above 
are H = 1, J q = 5(sj), l ql = S ql 5{t a )S(l - x/z) and 
S = S(ti — |f — q^)- F° r ne ld theoretic definitions of 
the beam, jet, and soft functions we refer the reader to 
[2] and references within. Large logarithms in the ratio 
of the scales hh ~ Pjt, ~ V PjtTi , and fis ~ t\ are 
resummed using the RG equations in SCET. 

In phenomenological studies, theoretical uncertainties 
to Eq.( 11 ) will arise from a truncation of the perturbative 
series in the hard, jet, and soft functions, the power cor- 
rections discussed earlier, and from non-perturbative soft 
function effects when t x < Aq C d- In addition, Eq.([TT]) 
will be affected by the standard PDF uncertainties. In 
the following, we make numerical predictions for the tj- 
distributions for the simplest case of a proton target. We 
estimate the perturbative uncertainties via a standard 
scale variation. We study the effects of non-perturbative 
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FIG. 1. da/ao = dydp^drx as a function n for a proton 
target at NLL accuracy. The bottom and top bands cor- 
respond to yfs = 90 GeV, P JT = 20 GeV, y = (Stage I 
EIC) and = 300 GeV, P JT = 20 GeV, y = (HERA) re- 
spectively. The dashed curve shows the singular part of the 
NLO cross-section for HERA kinematics. The scale choices 
for hh , Hj,^s are explained in the text. 
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FIG. 2. da/ao = dydp^drx as a function n for a proton 
target at NLL accuracy including non-perturbative Ti values 
with yfs = 300 GeV, P JT = 20 GeV, and y = 0. The different 
curves with peak positions from left to right correspond to 
(a, b, A) = (2.0, -0.2, 0.2), (1.2, -0.1, 0.3), (2.2, -0.4, 0.5), 
(1.8,-0.05,0.4) in Eq. |l5f respectively. 



soft radiation by modeling the soft function in the re- 
gion ti < A-qcd and varying the model parameters as 
described below. Numerical studies of power corrections 
are left for future work. Any deviations between data and 



the prediction of Eq. (11), after including the remaining 



uncertainties, will be a measure of the size of power cor- 
rections. These power corrections can be studied as a 
function of {t\,d, Pjt)- in the ratio in Eq. ([2|, the un- 
certainties from nuclear-independent effects will largely 
cancel out so that a deviation from unity probes nuclear- 
dependent effects. 

Fig. [I] shows the ri-distribution at Next-to-Leading- 
Log (NLL) accuracy for typical Stage-I EIC and HERA 
kinematics. We follow the conventions in Table 1 in 
Ref. [T5] for counting logs, set rif = 5, and use CTEQ6L 
PDFs [16 . The shaded bands correspond to NLL scale 
variation by choosing /i# = rPjT,H,j = "f\JP.]T T \-,^s = 
r t\ for the range r = {1/2,2}. The curve in the mid- 
dle of each band corresponds to r = 1. We also show 
(dashed curve) the singular part of the NLO cross-section 
for HERA kinematics. The effect of resummation is to 
tame the singular behavior of the fixed order cross-section 
in the T\ — > limit. The ri-distribution is cutoff at 1 GeV 
so that the soft function is still perturbative. 

In the region t\ ~ Aqcd, the soft function S{j\,\i) 
becomes non-perturbative and is modeled as 



= J dk a dkj5{Ti - k a - kj)S(k a ,kj,fi),(12) 



where S(k a ,kj 1 /i) is the generalized hemisphere soft 
function [2] . Non-perturbative physics in this hemisphere 
soft function is encoded by the convolution of the par- 
tonic soft function with a model function 5 mo d |17H19j 



as 



S(k a ,kj,fi) = / dk^dE S mod (k' a , k'j) 



xS pa . Tt .{k a - K,kj - k'j,^), (13) 

where <S par t. corresponds to the perturbative soft func- 
tion. S mo d(k' a ,k' b ) is chosen to peak near k' a b ~ Aqcd 
so that for t\ ^> Aqcd the soft function reduces to the 
perturbative result <S part .. For NLL resummation, the T\ 
distribution is sensitive to the non-perturbative physics 
only through the combination 



F mo d(u) 



S 



mod 



u + c 



c 



(14) 



This result is derived by expanding Eq. (11 1 to the NLL 
level and using u = k' a + hi, C = k' a — k' b in Eq. ( 13 ). 



The function Fmod is parameterized as 



fmod(«) 



(u - b) 2 

A 2 



,(15) 



where the normalization factor N(a, b, A) is chosen to 
satisfy J Q duF nm d(u) — 1, equivalent to the condition 
/ dk a J dkjS mod {k a ,kj) = 1. In Fig. [2] we plot the 1- 
jettiness cross section as a function of T\ for different pa- 
rameter choices of a, 6, and A as explained in the caption. 
We made the scale choices /j-h — P.j T , Mj = \]t\Pj t , and 
[i s = Ti^/l + (r™™/^) 2 with 7i min = 1 GeV. The dif- 
ferent models exhibit different behavior for t\ ~ Aqcd 
but converge to the perturbative result for larger values 
of T\ as required. The universality of the soft function 
allows its extraction from lepton-nucleon collisions to be 
then used in scattering off larger nuclei. 

A similar analysis can be performed for other nuclear 
targets. We leave further studies of nuclear phenomenol- 
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ogy, higher order resummation, and power corrections for 
future work. 

We thank Xiaohui Liu and Frank Petriello for use- 
ful discussions. This work was supported in part by 



by the U.S. Department of Energy under Contract 
No. DE-AC02-98CH10886 (ZK and JQ) and by the U.S. 
National Science Foundation under grant NSF-PHY- 
0705682 (SM). 



[1] I. W. Stewart, F. J. Tackmann and W. J. Waalewijn, 
Phys. Rev. Lett. 105, 092002 (2010). 

[2] T. T. Jouttenus, I. W. Stewart, F. J. Tackmann and 
W. J. Waalewijn, Phys. Rev. D 83, 114030 (2011). 

[3] X. -N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 
1480 (1992); R. Baier, Y. L. Dokshitzer, A. H. Mueller, 
S. Peigne and D. Schiff, Nucl. Phys. B 484, 265 (1997); 
B. G. Zakharov, JETP Lett. 65, 615 (1997); U. A. Wiede- 
mann, Nucl. Phys. B 588, 303 (2000); M. Gyulassy, 
P. Levai and I. Vitev, Nucl. Phys. B 594, 371 (2001); 
X. -F. Guo and X. -N. Wang, Phys. Rev. Lett. 85, 3591 
(2000); X. N. Wang and X. F. Guo, Nucl. Phys. A 696, 
788 (2001). 

[4] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. 
Lett. 91, 072301 (2003); J. Adams et al. [STAR Collabo- 
ration], Phys. Rev. Lett. 91, 172302 (2003); K. Aamodt 
et al. [ALICE Collaboration], Phys. Lett. B 696, 30 
(2011); S. Cha trchyan et al. [CMS Collaboration], 
|arXiv:1202.2554| [nucl-ex] . 

[5] I. Vitev, S. Wicks and B. -W. Zhang, JHEP 0811, 093 
(2008); I. Vitev and B. -W. Zhang, Phys. Rev. Lett. 104, 
132001 (2010). 

[6] J. Thal er and K Van Tilburg, JHEP 1202, 093 (2012) 

[arXiv:1108.2"70T1 [hep-ph] ] . 
[7] Z. -B. Kang, A. Metz, J. -W. Qiu and J. Zhou, Phys. 

Rev. D 84, 034046 (2011). 
[8] C. W. Bauer, S. Fleming, D. Pirjol and I. W. Stewart, 

Phys. Rev. D 63, 114020 (2001); C. W. Bauer, D. Pirjol 

and I. W. Stewart, Phys. Rev. D 65, 054022 (2002); 



[9] D. Boer, M. Diehl, R. Milner, R. Venugopalan, W. Vo- 
gelsang, D. Kaplan, H. Montgomery and S. Vigdor et al., 
|arXiv:1108. 17131 [nucl-th] . 

[10] M. Hirai, S. Kumano and T. -H. Nagai, Phys. Rev. C 76, 
065207 (2007); K. J. Eskola, H. Paukkunen and C. A. Sal- 
gado, JHEP 0904, 065 (2009 ); D. de Fl orian, R. Sassot, 
P. Zurita and M. Stratmann, |arXiv:1112.6324| [hep-ph], 

[11] I. W. Stewart, F. J. Tackmann and W. J. Waalewijn, 
Phys. Rev. D 81, 094035 (2010) [arXiv:0910.0467| [hep- 
ph]]. 

[12] M. Luo, J. -W. Qiu and G. F. Sterman, Phys. Rev. D 

50, 1951 (1994). 
[13] See, for example: F. Gelis, E. Iancu, J. Jalilian-Marian 

and R. Venugopalan, Ann. Rev. Nucl. Part. Sci. 60, 463 

(2010); and the references therein. 
[14] J. -W. Qiu and G. F. Sterman, Nucl. Phys. B 353, 137 

(1991). 

[15] C. F. Berger, C. Marcantonini, I. W. Stewart, F. J. Tack- 
mann and W. J. Waalewijn, JHEP 1104, 092 (2011) 

[16] P. M. Nadolsky, H. -L. Lai, Q. -H. Cao, J. Huston, 
J. Pumplin, D. Stump, W. -K. Tung and C. -P. Yuan, 
Phys. Rev. D 78, 013004 (2008) 

[17] Z. Ligeti, I. W. Stewart and F. J. Tackmann, Phys. Rev. 
D 78, 114014 (2008). 

[18] A. H. Hoang and I. W. Stewart, Phys. Lett. B 660, 483 
(2008). 

[19] S. Fleming, A. H. Hoang, S. Mantry and I. W. Stewart, 
Phys. Rev. D 77, 114003 (2008). 



